Nonlinear elastodynamics and energy transport in biopolymers 
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We present an exact analysis of the intrinsic dynamics of two types of biopolymers, described 
by the wormlike chain (WLC) and the wormlike rod chain (WLRC) models, respectively. We map 
each model to a classical field theory for an appropriate unit vector field defined on the space curve 
formed by the axis of the biopolymer. Conformation dynamics and energy transport are shown 
to be described, respectively, by the soliton-bearing Landau-Lifshitz equation and the nonlinear 
Schrodinger equation. This enables us to identify the soliton representing the "conformon" hypoth- 
esized in biopolymers. 
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In the years since the discovery of biological polymers 
such as proteins and DNA, it has been recognized [ij that 
the geometric shapes and topological properties of these 
long chain molecules play a significant role during pro- 
cesses such as replication and transcription. Experimen- 
tal techniques capable of studying the conformation and 
elastic properties of a single polymer by bending or twist- 
ing it are also available[2|. Semiflexible biopolymers such 
as actin have only a bending energy, and are described 
by the well-known wormlike chain (WLC) modelQ with 
a single elastic constant, the bending modulus. In con- 
trast, DNA with its double-helix structure is described 
by the wormlike rod chain (WLRC) model |^| with an 
additional elastic constant, the twist rigidity. Although 
considerable work has been done on various equilibrium 
properties of both these elastic models, their intrinsic 
dynamical properties have not been studied so far. The 
latter play a crucial role in the mechanisms of energy and 
information propagation along a biopolymer, an issue of 
vital interest to biologists, chemists and physicists alike. 

In this Letter we present, for the first time, a formal- 
ism to study the exact intrinsic dynamics of a polymer 
in the WLC and WLRC models. Several new results are 
obtained. By modeling the polymer as a slender elastic 
rod whose central axis is a space curve, the appropriate 
moving frame comprising the unit tangent, normal and 
binormal vectors is constructed for each model. We show 
that the WLC energy is a certain classical field theory 
Hamiltonian for the tangent vector, while the WLRC en- 
ergy is an analogous Hamiltonian for the normal vector 
field. Using a variational principle for the Hamiltonian, 
and imposing the unit vector constraint, the equilibrium 
conformations of the polymer are first found. Next, by 
setting the time derivative of the Hamiltonian to zero, the 
time evolutions of the vector fields consistent with their 
corresponding equilibrium conformations are found for 
the two models. It turns out that the equation of motion 
of the relevant vector field in each model resembles that 
of the precession of a magnetic moment 0, i.e., a "spin", 
in a self-induced "magnetic field". We also show that 



the dynamical equation in each of the models can be de- 
rived as Hamilton's equations from the respective elastic 
Hamiltonians, by using the z-component of the vector 
field and its azimuthal angle as the canonically conju- 
gate variables. As is well known, these variables give rise 
to the correct Poisson brackets for the "spin" compo- 
nents. Remarkably enough, the time evolution equations 
for the tangent vector in the WLC model and the nor- 
mal vector in the WLRC model have the same form as 
the well-known Landau-Lifshitz equation. The latter is 
completely integrable by the inverse scattering transform 
method. We show that the writhe and the twist of the 
vector field in each model are conserved quantities under 
time evolution. By White's theorem^], this leads to the 
conservation of its linking number as well. 

The non-extensibility of the polymer during time evo- 
lution entails compatibility conditions on the frame fields 
of the polymer. Imposing these leads to the completely 
integrable nonlinear Schrodinger equation (NLS) Q for 
the energy density in each model. This shows that en- 
ergy can propagate along the polymer like a localized 
wave packet, in both types of biopolymers. 

Using the one-soliton solution of the NLS equation, we 
are able to identify the curvature and torsion functions 
in the two models. From these expressions we find that, 
in the WLC model, the soliton corresponds to a local- 
ized "loop" propagating along the semiflexible polymer, 
while in the WLRC model (applicable to DNA), it is 
possible to have a localized "supercoiled" entity propa- 
gating along the polymer, for appropriate values of the 
width and velocity parameters of the soliton. We in- 
terpret our solutions as the explicit manifestation of the 
so-called "conformon", an entity that has been hypothe- 
sized by various authors 0, 0, [HI to play an important 
role in biology - for example, these are localized packets 
of conformational strain on mitochondria. 

In the context of an cv-helical protein, the NLS soliton 
first found by Davydov|ll} using a quantum mechanical 
lattice model is well studied. However, it has been 
noted^3| that the Davydov soliton is essentially a 
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classical object. Our analysis leads directly to a classical 
soliton. It also enables us to suggest a mechanism for 
charge transport along the polymer, in the light of the 
results of Goldstone and Jaffe 13J on the curvature- 
induced localization of a quantum particle propagating 
in a thin curved tube. 

The WLC model: A semifiexible biopolymer is repre- 
sented as a slender elastic rod with a single elastic mod- 
ulus characterizing the bending energy of the polymer 
. In general, the axis of the polymer is a space curve. 
The unit tangent t, normal n and binormal b = t x n 
comprise a frame on the space curve. They satisfy the 
well known Frenet-Serret equations^) t s = sn , n s = 
— Kt + rb , b s = — m. Here s is the arclength and the 
subscript denotes the corresponding derivative. As usual, 
k is the curvature and r the torsion of the curve. The 
Frenet-Serret equations can be written in the compact 
form F s = O x F where F stands for t, n or b, and 
f2 = t t + k b. Thus, as one moves along the curve, the 
frame rotates with a "spatial angular velocity" fl. 

Consider a polymer with a constant disc-shaped cross- 
section. We may imagine the polymer to be formed 
by the tumbling of this disc through space ^| with- 
out twisting about its axis, thereby implying that its 
angular velocity should always lie in the plane of the 
disc. This amounts to the Fermi- Walker parallel trans- 
port of the triad, as has been recognized in another 
context [l5|. In order to satisfy this condition, we in- 
troduce new unit vectors N = n cos a(s) + b sin a(s) 
and B = — n sin a(s) + b cos a(s) , where a(s) — 
— J_ r(s') ds' . Then the triad comprising (t, N, B) can 
be shown to satisfy the Fermi- Walker (FW) equations 
t s = k (N cos a — B sin a) , N s = — t k cos a , and 
B s = t k sin a . These equations lead to an angular 
velocity u) — k(N sin a + B cos a) for the triad. The 
corresponding Hamiltonian is therefore of the form 



H = 



A 



lu ds 



k 2 ds 



t c ds, 



(1) 



where A is a constant. This is precisely the expression 
written down for the elastic energy of the rod in the WLC 
model, with A standing for the bending modulus. 

We first obtain the minimum energy conformations of 
the WLC. We scale out A and use the variational princi- 
ple SH = on the Hamiltonian 
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(3(t 2 - l)}ds 
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where (3 is the Lagrange multiplier corresponding to the 
unit vector field constraint t 2 = 1. A short calculation 
yields 



t x t.„ 



0. 



(3) 



whose solutions represent the static conformations 
sought. To find the dynamics of the polymer, we set 



t = t(s.u) in Eq. Q where u denotes the time. The 
energy conservation condition dH/du = then leads to 
J t su • t s ds — , which yields, on integrating by parts, 
J t 8S ■ t u ds = 0. In turn, this gives t u = V x t ss , where 

V is an arbitrary vector. The static solution of this equa- 
tion is V x t ss = 0. A comparison of this with Eq.JHJl im- 
mediately suggests the simplest possible choice, namely, 

V = t. The corresponding time evolution for t is then 
given by 

t tt = t X t ss . (4) 

This equation resembles that 5] for the precession of 
a magnetic moment or spin in a self-induced mag- 
netic field. To make this connection more precise, 
we must ask whether Eq. Q can be obtained di- 
rectly from the WLC elastic Hamiltonian in Eq. 0J. 
Now, the unit vector t may be parametrized as t = 
(sin 9 cos 4>, sin 9 sin cj>, cos 9). If t is to be a classical 
spin vector, its cartesian components must satisfy the 
standard Poisson bracket relations {t a , t@} = e OL ^ 1 t 1 . (In 
turn, this is ensured if we have {0, cos 9} — 1.) With 
this input for the canonical Poisson brackets, we can show 
that Eq. can be obtained from the equation of motion 
t" = {t a ,H} corresponding to the WLC elastic Hamil- 
tonian. Further, the canonical Poisson bracket relations 
also enable us to show that the writhe of t, defined by 
Wr = (2tt)~ 1 J <j) s (1 — cos 9)ds , is a conserved quan- 
tity, i.e., {H, Wr} — 0. The writhe is a measure of the 
chirality of the curve. 

We note that Eq. (@J, which is called the Landau- 
Lifshitz (LL) equation in the spin context, is satisfied by 
the tangent t along the axis of the polymer in the WLC 
model. The LL equation appears as a spin evolution 
equation in the dynamics of a continuum Heisenberg fer- 
romagnetic chain. It is a completely integrable equation 
with strict soliton solutions and an infinite number of 
constants of motion. Next, we use the FW equations for 
(t,N,B) in Eq. (0J, and impose the compatibility con- 
ditions t su — t us , B stI = B us . After some algebra, we 



get k u 
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k, r and r„ 
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These equations are just the imaginary and real parts 
of the nonlinear Schrodinger equation (NLS) 



ts[llj 



o. 



(5) 



when ifj = k exp (i J r ds) is the well-known Hasimoto 
function [l^. (In the present context, the compatibil- 
ity conditions represent the physical property of non- 
extensibility of the polymer.) Thus | - 0| 2 = k 2 , the energy 
density of the polymer. From the equation for t u , we can 
easily show that the polymer evolves in such a way as to 
conserve its total twist, namely, Tw — J r ds. Since we 
have already shown that the writhe Wr is conserved, it 
follows from White's theorem Tw + Wr = Lk, that 
the linking number Lk of the WLC polymer is also con- 
served. 
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The one-soliton solution of Eq. (JSJ) is given by 

ip(s, u) = pi sech^i exp (i£ 2 ) , (6) 

with £i = pi(s - v x u) and £ 2 = \v\s - \(v 2 - p\)u, 
where the amplitude p\ and the velocity v\ of the en- 
velope of the solution are independent parameters. The 
energy density is | ip | 2 = p\ sech 2 £i . This represents an 
energy lump propagating with velocity v\ along the poly- 
mer. For this solution the curvature n = p\ sech£i and 
the torsion r = It is clear from these expressions 

that this solitonic excitation represents a localized defor- 
mation of the curvature, with localization length p^ 1 , 
that moves along a helical polymer with a constant pitch 
\v\. We believe that the foregoing provides a precise and 
quantitative dynamical underpinning for the conformon 
concept put forward by Green and Ji and by others 0. 

We can go even further. Let us first consider the spe- 
cial case of a planar polymer, by setting r = in the FW 
equations. We get t s = kN , N s = — nt , and B = B c , 
a constant vector. Then, using the LL equation, a bit 
of algebra gives t u = ft s B c and N„ = (k ss /k)B c . The 
non-extensibility condition t su — t us yields k u = , 
implying that n = n[s). The compatibility condition 
N su = N us gives, after the rescaling s — > 2~ 1//2 s, the 
equation k ss + \k 3 = Ek where E is an arbitrary 
constant. But this equation is just the time-independent 
counterpart of the time-dependent Schrodinger equation 
itpu + V'ss — Vip = for a particle with wave function 
ip = k exp(iEu), in the "potential" V = — jK 2 (s). In- 
terestingly, it has been shown by Goldstone and Jaffe[l3j 
that the one-dimensional Schrodinger equation for a 
free particle in a thin tube whose axis follows a space 
curve of curvature k acquires an effective attractive 
potential -|k 2 (s) (for small k). In our formalism, 
this has emerged naturally as the elastic energy of 
the WLC. The time- independent Schrodinger equation 
above has a (static) soliton solution k — p sech ps 
where p = E 1 / 2 . Combining these results, we have 

V = — jp 2 sech 2 ps , and a corresponding localized wave 
function ip — p sech (ps) exp(iEu). Applying this result 
to an electronic wave function on a polymer, this suggests 
the possibility that a soliton, i.e., a nonlinear excitation 
of the elastic polymer, can lead to electron localization. 
In the non- planar case (r ^ 0), the effective potential 
continues to be V = — \k 2 , but now k — k(s,u). In 
particular, we can now interpret the wave function ip as 
the one-soliton solution in the time-dependent potential 

V = — \p\ sech 2 £i. In the context of electron localiza- 
tion, the electron gets trapped in a moving potential 
well that travels along the polymer. This ties in with 
the conformon concept as su gge sted by VolkensteinQ 
and by Kemeny and Goklany[lfJ. We conjecture that 
this moving soliton plays a role in the unexpectedly high 
electrical conductivity found in biological membranes. 



The WLRC model: Here a polymer such as DNA is rep- 
resented as a slender twisted elastic rod with two elas- 
tic constants that characterize the bending energy and 
the twisting energy, respectively^. The state of the 
polymer is specified by the configuration of the ribbon 
formed by a longitudinal cross-section of the rod con- 
taining its central axis. The ribbon has a uniform twist 
angle luq per unit length around the central axis, to take 
into account the double helical structure of DNA. In gen- 
eral, the central axis of the ribbon is a space curve. 
Let t be the unit tangent to this axis at any point, 
and No the unit normal to the surface of the ribbon 
at that point. Defining B = t x N , the axis be- 
comes a framed space curve on the surface of the ribbon. 
It is obviously convenient to work with a "co-twisted" 
frame (t,N,B) where N = N cos ujqs + B sin u s 
and B = —No sin ojqs + Bo cos luqs . In this frame the 
axis of the polymer is a space curve on the surface 
of the ribbon, and hence is described by the Darboux- 
Ribaucour (DR) equations [T^ t s = kn N + K g B , N s = 
— kn t + T g B , B s = —Kg t — Tg N. Here kn and K g rep- 
resent, respectively, the normal and the geodesic curva- 
tures, and T g is the geodesic torsion. As the width of the 
ribbon goes to zero, the geodesic curvature of any contin- 
uous curve on its surface tends to zero. The only curve 
which survives is the axis of the ribbon, with curvature 
kn and torsion r g . The DR equations reduce to 

t s = k w N , N s = -K N t + T g B , and B s = -r ff N , (7) 

where N and B now denote the normal and binormal to 
the curve. Once again, as described earlier in the case 
of the WLC model, the frame (t,N,B) can therefore be 
regarded as rotating with a "spatial angular velocity" 
$7 = T g t + kjvB. Hence the corresponding Hamiltonian 
is of the form 

H= l - J(AK 2 N + CT 2 )ds=^ J(At s 2 + CB s 2 )ds, (8) 

where A and C are constants. But this expression is 
precisely the elastic energy of a twisted polymer in the 
WLRC modelHEJ, where A is the bending modulus as 
before, and C is the twist modulus. 

We first discuss the symmetric case A = C. Using Eq. 
the (scaled) energy can be written as H = | J N 2 ds. 
We proceed as in the case of the WLC model. The con- 
straint in this case is N 2 = 1. We finally obtain 

N u = N x N ss (9) 

as the counterpart of Eq. Q. The discussion following 
the latter equation applies to the present case, with N 
replacing t. Thus it is now N that behaves like a spin 
vector field. 

We drop the subscripts in kn and r g hereafter, to sim- 
plify the notation. Equations (7J and JjJJ give N u = r s t+ 
k s B. This, along with the conditions N-t = , N-B = 
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and t su = t us , yields t u = — r s N + ^(k 2 + r 2 )B and 
B M = —\{k 2 + r 2 )t — k s N. The compatibility condition 
N s „ = N us then leads to[20| 

K„+T s;j + ^(K 2 +r 2 ) =0, r tt -K ss -|(/« 2 +r 2 ) =0. (10) 

A short calculation shows that Eqs. (|10|l represent re- 
spectively, the real and imaginary parts of the NLS, Eq. 
|[SJ|, if we set ip = k + it. In the WLC model, in contrast, 
it is the function k exp [i J r ds) that satisfies the NLS. 

For the one-soliton solution of the NLS in Eq. ©, the 
energy density in WLRC model once again has the form 
\ip | 2 = p 2 sech 2 £;i , representing an energy lump propa- 
gating with velocity v\. However, the curvature and tor- 
sion differ from those in the WLC model. They are now 
given by k — pisech£i cos £2 and r = pisech£i sin £2- 
In contrast to the WLC model, the curvature and torsion 
are both localized functions, with an additional spatial 
periodicity A = Anv^ 1 . This is the signature of a prop- 
agating localized "supercoiled" excitation (a conformon) 
in DNA. 

Finally, we consider the general case A ^ C, for which 
the scaled Hamiltonian is H = i j(n 2 + f 2 ) ds where 
f = Ct/A. The NLS is now satisfied by the function ip = 
k + if. The results for the one-soliton solution described 
in the preceding paragraph continue to apply. The only 
difference is that the torsion of the polymer is now scaled 
by a factor A/C. 

To conclude, we have provided an exact analysis of 
the nonlinear dynamics of elastic polymers in the WLC 
and WLRC models using a geometric formulation, and 
obtained several new results. Although we have consid- 
ered only open polymers, our analysis can be extended 
to closed polymers as well. As a consequence of the pe- 
riodic boundary conditions in this case, the solutions of 
the NLS now turn out to be elliptic functions. 

Based on our results we note the following similarities 
between our soliton and the "conformon" conjectured in 
molecular biology 0, 0: (i) It is an energy packet as- 
sociated with a conformational strain localized in a re- 
gion much shorter than the length of the molecule; (ii) 
it needs a "rigid framework" (i.e., a structure with finite 
bend/twist moduli) for its formation; and (hi) its prop- 
erties appear to depend on the coupling between adja- 
cent bonds, e.g., in a polypeptide chain. The compati- 
bility of our results with (i) and (ii) is obvious. As for 
(iii), we note that in both the WLC and WLRC mod- 
els, the Hamiltonian density involves spatial derivatives 
of unit vector fields of the form w 2 . But this is just 
the continuum limit of the nearest-neighbor interaction 
w, • Wi + i where i is the site index on a discrete chain. 
Further, as we have seen, w(s) acts like a spin vector 
field. Thus the dynamical evolution of the polymer is 
governed by the nearest-neighbour (i.e., adjacent) cou- 
pling of an effective "spin vector" at each site on the 
chain. As is well known, there are several mechanisms in 



condensed matter physics which can lead to such effective 
spin Hamiltonians. A secondary effect can also occur, as 
we have discussed. As our soliton (localized curvature) 
moves on the polymer, it can trap an electron and drag 
it along. This is in accordance with the conformon con- 
cept of VolkensteinQ and of Kemeny and Goklany[T3. 
This mechanism may be important in the understanding 
of charge transport in biological molecules, where cur- 
vatures abound. As we have demonstrated, curved ge- 
ometry and nonlincarity are intimately related, and we 
believe that our results are also of significance in other 
kinds of transport phenomena in molecular biology. 
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